Phase-ordering kinetics: ageing and local 
scale-invariance 



Malte Henkel 

Laboratoire de Physique des Materiaux ( CNRS UMR 7556), Universite Henri Poincare Nancy I, 
B.P. 239, F - 54506 Vandceuvre-les-Nancy Cedex, France 

Abstract. Dynamical scaling in ageing systems, notably in phase-ordering kinetics, is well- 
established. New evidence in favour of Galilei-invariance in phase-ordering kinetics is described. 



DYNAMICAL SCALING IN AGEING SYSTEMS 

The study of the long-time dynamics of statistical systems far from equilibrium has been 
a topic of intensive study. In many instances, the relaxation times towards equilibrium 
can become extremely long such that the system stays for all intents and purposes out of 
equilibrium. A paradigmatic example is provided by glassy systems which might be con- 
sidered as an extremely viscous liquid. In principle, the presence of very long relaxation 
time-scales might suggest that quantitative properties of glassy dynamics should depend 
on a huge variety of microscopic 'details' and furthermore, as the behaviour of a glass 
may depend on its previous (thermal, mechanic,. . . ) history, those system age. However, 
as first pointed out by Struik in 1978 [1], the ageing dynamics of many physically very 
different glass-forming systems can be described in terms of universal master curves. 
The challenge is to try to understand the origin of this dynamical scale-invariance and to 
compute the form of these master curves from the essential characteristics of the system. 

Here we shall consider an analogous situation in supposedly more simple systems 
without intrinsic disorder or frustration. From some initial state (typically one chooses a 
fully disordered initial state) the system is rapidly 'quenched' either onto its critical 
point or else into a region of the phase diagram with at least two stable stationary 
states. The dynamics and an eventual ageing behaviour (that is, a breaking of time- 
translation invariance) is then observed. One distinguishes physical ageing, where the 
underlying microscopic processes are reversible, from chemico-biological ageing, where 
irreversible microscopic processes may occur. In terms of models, a simple example for 
physical ageing is provided by the phase-ordering kinetics of a simple Ising ferromagnet 
with purely relaxational dynamics quenched to below its critical temperature T c > 0. On 
the other hand, the ageing behaviour of the contact process (particles of a single species 
A move diffusively on a lattice and react according to A +A — ► and A — ► 2A) provides 
a paradigmatic case of ageing with underlying irreversible processes. 

Physically, these two kinds of ageing phenomena are quite distinct. In relaxing fer- 
romagnets, see e.g. B2L ]§\ for reviews, there is for T < T c a non-vanishing surface 
tension between the ordered domains which leads to the formation and growth of or- 



dered clusters of linear size L — L{t) ~ t x l z , where z is the dynamical exponent. For 
purely relaxational dynamics, it can be shown that z = 2 for T <T C whereas for T = T c , 
the non-trivial value of z equals the one found for equilibrium critical dynamics. If (t, r) 
denotes the time- and space-dependent order parameter, it is convenient to characterize 
the ageing behaviour in terms of the two-time correlation and linear response functions 

(1) 

where h is the magnetic field conjugate to 0. The autocorrelation and linear autoresponse 
functions are given by C(t,s) = C(t,s;0) and R(t,s) = R(t,s;0). Here and later space- 
translation invariance will be assumed. In phase-ordering, dynamical scaling is found in 
the regime where 

t 3> T m i cro , s 3> T m i cro and t s 3> T m i cro (2) 

where Tmicro is some microscopic reference time-scale. In the ageing regime © the only 
relevant lengths scales are describes in terms of L(t) and one expects 

C(t,s)=s- b f c (t/s) , R(t,s) = s - l - a f R {t/s) (3) 

with the asymptotics fc,R{y) ~ y - ^/ 2 for y large. We stress the importance of the third 
condition in eq. © for the validity of the scaling forms ©. Throughout, the scaling limit 
t,s — > °° with y = t/s> 1 fixed will be implied. While the autocorrelation (autoresponse) 
exponents Xc,r are new, independent exponents, the exponents a,b can be explicitly 
given. At criticality T = T c , one has a = b= (d — 2 + r\)/z = 2/3/ vz, where /3, v,r/ are 
standard equilibrium critical exponents. For T < T c , one has b = always. In simple 
scalar systems with short-ranged equilibrium correlations, such as the Ising or Potts 
models (with d > 1), one has a = 1/z = 1/2. In the case of long-ranged equilibrium 
correlations, a may be different, i.e. a = id — 2) jz in the J-dimensional spherical model. 

On the other hand, in irreversible ageing systems such as the critical voter-model 
or the critical contact-process, there is no surface tension and the dynamics proceeds 
through cluster dissolution ||5|]. Still, in the ageing regime © one observes again the 
same formal scaling behaviour eq. ©. However, the exponents a and b need no longer 
be the same even if one considers the ageing at a phase-transition in the non-equilibrium 
steady-state. For example, there is recent numerical evidence from the critical contact 
process in both ID and 2D that b = 2f3/v±z which naturally generalizes the result of 
critical systems with detailed balance, but a = b — 1 



C(t : s;r) = (0(/\F)0(s,O)) , R(t,s;r 



Sh(s,0) 



GALILEI-INVARIANCE IN PHASE-ORDERING KINETICS 

Having reviewed current knowledge on the dynamical scaling of ageing systems, no- 
tably on the values of the ageing exponents a and b, we now consider the scaling func- 
tions fc,R(y) themselves. In particular, we ask whether there exist any generic, model- 
independent argument which might inform us about the form of the functions fcjt(y)- 



In this context, the following general statements about the dynamical symmetries of 
phase- ordering kinetics can be made. 

1 . Time-translation invariance is broken. 

2. There is dynamical scaling an, i.e. formally the order parameter satisfies the 
covariance condition 

(j)(t,r) = a x ^(j)(a 2 t,ar) (4) 

where a is a constant rescaling factor and x§ an exponent. 

3. There is new evidence for Galilei-invariance which we now discuss. 

Conventionally, one starts from a coarse-grained order parameter which is assumed 
to satisfy a Langevin equation, for example the one for model A dynamics J2] 

2^f = A0-— -^ + 7] (5) 
dt d0 

where T = (2^)~ l is a kinetic coefficient, A is the spatial Laplacian, V((j>) is a typical 
double-well potential (e.g. V (0) = ((j) 2 — l) 2 ) and r\ is the thermal gaussian noise with 
covariance (r\(t)r\(t')) = 2T8(t — t'). In addition, one assumes a gaussian uncorrelated 
initial state with covariance (0(0, r) 0(0,0)) = ao8(r). The associated field-theoretic 
action in the Martin-Siggia-Rose (MSR) formalism is 

S[0,0]=S o [0,0]+S,[0] (6) 

where is the response field and 



So [0,0] = J dtdr 



<$>(2J£d t - A) <$> + <$>- 
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S,[0] = -r|d?dr0 2 (?,r)-^|dr0 2 (O,r) (7) 
Then autocorrelators and autoresponses can be found as follows 

C(t,s) = <0(O0Os)> , R(t,s) = U{t)Hs)) (8) 



We shall consider as an extension of dynamical scaling the transformations of the 
so-called Schrodinger group, defined by ISD 

at + B ^ &r + vt + a 

t -> F > r ~* 5 ( 9 ) 

yt + 8 yt + 8 

where a8 — j3y = 1 and M is a J-dimensional rotation. Besides translations in time 
and space, rotations and scale transformations with z = 2, the Schrodinger group also 
contains the so-called 'special' transformations parametrized by y. The Schrodinger 
group is the maximal dynamical symmetry group of the free (and also of several non- 
linear) Schrodinger/diffusion equations and acts projectively (i.e. up to a phase factor) 
on the wave function 



The relationship of the Schrodinger group to phase-ordering kinetics can now be 
formulated in terms of the following three theorems. 

We consider an arbitrary space-time infinitesimal coordinate transformation dr^ = 
with /i = 0, 1, . . . ,d. Let rj stand for the phase picked up by the wave function under 
such a transformation. We call a MSR-theory local if the MSR-action © transforms as 

8S= [dtdr (r Mv a M e v +/ M a M i])+ / dr (U v e v +VJj) (10) 
J ' ' J(t=0) 

Here T^v is the energy-momentum tensor, a conserved current and the second integral 

is only over the initial line t = 0. 

Theorem 1: [10] For a local MSR-theory, one has 

phase-shift invariance 
space-translation invariance 
scale -invariance with z = 2 
Galilei-invariance 

This result is completely analogous to a well-known result in conformal field-theory. We 
point out that the requirement of time-translation invariance is not required in order to 
derive the invariance under special Schrodinger transformations for local theories. 

To prove this, it suffices to write down the various infinitesimal transformations ex- 
plicitly. From invariance under phase shifts it follows V = 0, space-translation invariance 
implies U\ = . . . = Uj = 0, from scale-invariance it follows 2Tqq + T\ \ + . . . + = and 
finally Galilei-invariance implies 7$ + 2^Ji = for i = 1 , . . . , d. Consequently 8S = 
under special Schrodinger transformations, as asserted. q.e.d. 

One may write down the tensor T^ v and the current explicitly, for free fields this 
has been done in ifioll . 

We now consider the role of Galilei-invariance more closely. It is well-known that 
a system coupled to a heat bath with a uniform temperature T > cannot be Galilei- 
invariant. This can be easily seen from the MSR-action ©, since the noise terms Sb[<j>] 
are not invariant under a phase-shift — >• — » e _r, 0. At most, the deterministic 
part So of the action might be Galilei-invariant. If that is the case, one has the following 
Bargman superselection rule 

( ) ~ 8, hm (12) 

n m 

The index refers to the average taken with respect to the deterministic part Sq only. 
Theorem 2: [11] If the deterministic part Sq of the MSR-action (0) 15 Galilei-invariant 
such that A12\l holds true, then 



R(t,s) = 
C(t,s) = 




where the two-point function R Q (t,s) = (^)(t,y)^)(s,y))o and the three-point function 
R (t,s,u;r) = ((j>(t, y)(j>(s, j)0 2 (m,F-|-j))o are fixed by the deterministic part Sq. 




special Schrodinger invariance (11) 



To see this, one merely has to include the noise term S&[0] into the average, viz. 

R = (00) = (00e~ 5fc W)o = (00)o because of the Bargman superselection rule (IT2b . 
C(t, s) is found similarly. q.e.d. 

Consequently, given the Galilei-invariance of the deterministic part, one can find the 
contributions of both the thermal and the initial noise. Remarkably, the response function 
is noise-independent, whereas the correlator vanishes in the absence of noise. 

Finally, we have to address the question whether the deterministic part of the Langevin 
equation © can be Galilei-invariant. At first sight, the answer seems to be negative, 
since a well-known mathematical fact [9] states that the non- linear Schrodinger equation 

(2mid,-A)0=F(*,r,0,0*) (15) 

is Schrodinger-invariant only for the special potential F = c(00*) 2 / d 0, where c is a 
constant. Furthermore, the solutions are necessarily complex. These difficulties can 
be circumvented by considering the mass m as a new dynamical variable lIlzllOll . We 
introduce a new wave function i// via 

/CO 
d;e- im ^ ¥ (C,t,r) (16) 
-00 

and look for the symmetries of the new non-linear 'Schrodinger' -equation 

(2d^ t -d^)xif = gF{^t^xf) (17) 

where we wrote the coupling constant g of the non-linear term explicitly. Since from 
dimensional counting, g is in general dimensionful, it should transform under the action 
of the Schrodinger group. We have systematically constructed all such representations of 
the Schrodinger Lie algebra and also of its subalgebra when time-translations are left out 
JrHl . We then find all invariant non-linear equations of the type (flTt . The full analysis 
will be presented elsewhere, here we merely quote one special result. 
Theorem 3: E^JIfwe writethe solutions of ^3 W= Vg(C,t,r) = g^-^l^^^.t.r), 
where x is the scaling dimension of y/ and y the scaling dimension of g, then the equation 

{2d c d t - d}) m = g-^yf (V / 4 V) (i8) 

is Schrodinger-invariant, where f is an arbitrary function. 

We have obtained in this way a formulation of Schrodinger-invariance which allows 
for real-valued functions and furthermore should be flexible enough to include the 
double- well potentials which enter into the Langevin equation ©. 

TESTS OF GALILEI-INVARIANCE 

The results quoted above in the theorems 1-3 lead to quantitative predictions which have 
been successfully tested in simulations. We first consider the response function, which 




FIGURE 1. Scaling function p^\y,jx) in the 2D Ising model at T = 0.66 T c for two values of the control 
parameter fi, as a function of y = t /s (after IU4I0 . 



according to theorem 2 can be found from the assumption of its covariant transformation 
under the Schrodinger group. This leads to 

/ M f 2 \ 

R(t : s;r) = R( t , s )expl- — — J , R(t,s) = s^Ut/s) 

f R (y) = foy 1+a '- x " /z (y-ir 1 - a ' (19) 

where <M and /q are non-universal constants. In most cases, one has a = a', but ex- 
ceptions are known to occur, e.g. in the ID Glauber-Ising model. The gaussian space- 
dependence of R(t, s;~r) is characteristic of Galilei-invariance. Since response functions 
are very much affected by noise and hence difficult to measure directly, a convenient 
way of testing (IT9l) is to study the spatio-temporally integrated response 

[°du I ^drT 4 - 1 R(t,u;r)=r s d/2 - a p ( - 2 \t/s,n) + --- (20) 
Jo Jo 

where /i is a control parameter, ro a constant related to fo and the function p( 2 ) can be 
found explicitly from dT9t [ 14J]. 

In figure [T] we compare the prediction derived from (TH21) with simulational data 
ifLHl in the 2D kinetic Ising model, quenched into its ordered phase and with a purely 
relaxational heat-bath dynamics. We used the exponents z = 2, a = a' = l/z = 1/2 and 
Xr = 1.26. From the plots, one sees a nice collapse of the data obtained for several 
values of the waiting time s. Finally, the full curve agrees perfectly with the data. Similar 
results also hold true in the 3D Ising model llbUl and also for the 2D three-state Potts 
model lfl5ll . Furthermore, the prediction ( fT^l) can be reproduced in the exactly solvable 
J-dimensional spherical and ID Glauber-Ising models. 

We remark that R(t, s) is well reproduced in the critical ID contact-process JfJ. 



As a second example, we consider the calculation of the autocorrelator C(t,s). From 
theorem 2, this requires the calculation of a noiseless three-point response function. 
Since Schrodinger-invariance fixes the three-point function only up to an undetermined 
scaling function, a further argument is needed in order to determine C(t,s) completely. 
For models which are described by an underlying free-field theory, it turns out that 
a simple heuristic idea based on the absence of singularities in C(t,s) leads to the 
following simple expression for T = in the scaling limit ifTlll 

C(t, S )=f c (t/s) , fc(y)^C ((y+\) 2 / y y Xc/2 (21) 

This agrees indeed with the exact solutions of the spherical model, the spin-wave 
approximation of the XY model, the critical voter model and the free random walk. 
It can also be checked 111 ill that the terms coming from thermal noise are irrelevant, as 
expected from renormalization group arguments yfl. 

Finally, for models not described by a free-field theory one can invoke an extension 
of Schrodinger invariance to a new form of conformal invariance E.ThenCXM) can 
be written in terms of hypergeometric functions and this prediction again agrees nicely 
with simulational data in both the 2D Ising IH6ll and three- states Potts models 11511 . 

Summarizing, we have presented evidence, both conceptual and simulational, which 
strongly suggest that dynamical scaling in phase-ordering kinetics can be extended to 
the larger Schrodinger group (without time-translations) of local scale-transformations. 
In particular, it appears that equations of the form (fl"8l) . rather than the Langevin equa- 
tion © used traditionally, allow for a simple symmetry characterization. However, the 
derivation of equations of the type (IT~8l> from a physical argument remains to be under- 
stood. 
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